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A novel approach is proposed for realizing an arbitrary- 
operation gate with a SQUID qubit via pulsed-microwave ma- 
nipulation. In this approach, the two logical states of the 
qubit are represented by the two lowest levels of the SQUID 
and an intermediate level is utilized for the gate manipulation. 
The method does not involve population in the intermediate 
level or tunneling between the two logical qubit states during 
the gate operation. Morever, we show that the gate can be 
much faster than the conventional two-level gate. In addition, 
to take the advantage of geometric quantum computing, we 
further show how the method can be extended to implement 
an arbitrary quantum logic operation in a SQUID qubit via 
geometric manipulation. 

PACS numbers: 03.67.Lx, 03.65.-w, 74.50.+r, 85.25.Dq 



In the growing field of quantum computing, supercon- 
ducting devices including cooper pair boxes, Josephson 
junctions, and superconducting quantum interference de- 
vices (SQUIDs) [1-10] have appeared to be among the 
most promising candidates for quantum computing. Rea- 
sons for this are that these systems can be easily fab- 
ricated to implement a large-scale quantum computing 
and have been demonstrated to have relatively long de- 
coherence time [11,12]. In the past few years, for SQUID 
systems, many methods for demonstrating macroscopic 
quantum coherence [7], creating a superposition state 
(i.e., a Shordinger cat state) [9], or completing pop- 
ulation transfer between low-lying levels (i.e., a NOT 
gate) [10] have been presented. However, the key point 
for quantum computing is how to realize an arbitrary 
single-qubit operation plus a two-qubit controlled-NOT 
(or controlled-phase shift) operation, since this set of op- 
erations makes a universal quantum computing [13]. 

In this letter, we focus on how to achieve an arbitrary- 
operation gate with a SQUID qubit. The proposed 
method has several advantages: (a) during the gate op- 
eration, the intermediate level is not populated and thus 
decoherence induced by spontaneous emission from the 
intermediate level is greatly suppressed; (b) no tunneling 
between the SQUID-qubit levels |0) and |1) is required 
during the gate operation, therefore the decay from the 
qubit levels can be made negligibly small via increasing 
the potential barrier between the qubit levels; (c) more 



importantly, as shown below, the gate operation can be 
performed faster compared with the conventional gate 
operation. In addition, we show that the approach can be 
extended to perform an arbitrary operation on a SQUID 
qubit via geometric manipulation of phase. 

Consider a SQUID driven by two classical microwave 
pulses / and II. In the following, the SQUID is treated 
quantum mechanically, while the microwave pulses are 
treated classically. The Hamiltonian H for the coupled 
system can be written as 



H = + H 



H 



(ii) 
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where H s is the Hamiltonian for the SQUID; and H\ 
and ff} are the interaction energies for the SQUID 
with the two pulses, respectively. 

The SQUIDs considered throughout this letter are rf 
SQUIDs each consisting of a Josephson tunnel junction 
enclosed by a superconducting loop. The Hamiltonian 
for an rf SQUID (with junction capacitance C and loop 
inductance L) can be written in the usual form 
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where <E> is the magnetic flux threading the ring, Q is 
the total charge on the capacitor, & x is the static (or 
quasistatic) external flux applied to the ring, and Ej 
= 7 c $o/27T is the Josephson coupling energy (I c is the 
critical current of the junction and $0 — h/2e is the 
flux quantum). In practice, the single junction is often 
replaced by a dc SQUID with low inductance I <C L which 
is effectively a junction with its critical current tunable 
by varying the flux applied to the dc SQUID. 

The SQUID is coupled inductively to the fcth pulse 
(k = I, II) through the interaction energy Hj* = 
\k ($ — $x) Q^w: where = — 1/L is the coupling co- 
efficient between the SQUID and the fcth pulse; and 
is the magnetic flux generated by the magnetic 
component B fe (i) of the fcth pulse, which is given by 
= J s B k (t) ■ dS (S is any surface that is bounded 
by the ring). The expression of B *-'(£) takes the form 
B fe (t) = Bgcos (ijJkt + 4>k) , where Bg, u>k and <f>k are the 
amplitude, frequency and phase of the fcth pulse, respec- 
tively. 

Consider a A-type configuration formed by the two 
lowest levels and an intermediate level of the SQUID, 
denoted by |0),|1) and \a) with energy eigenvalues 
E , Ei, and E a , respectively [Fig. 1(a)]. Suppose that the 
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coupling of |0) , |1) and \a) with the other levels via the 
pulses is negligible, which can always be realized by ad- 
justing pulse frequencies, or level spacings of the SQUID. 
Under this consideration, it is easy to find that (i) if no 
|0) <-» |1) and |1) <-» | a) transitions are induced by the 
first pulse, i.e., the first pulse is far-off resonant with the 
|0) «-> |1) and |1) <-» | a) transitions, which can be real- 
ized by setting |^oi — u>i\ 3> figi an< ^ ~ w /| ^ ^io> 
(ii) if the second pulse is far-off resonant with |0) «-> |1) 
and |0) «-> | a) transitions, i.e., |o»oi — oju\ ^> il^i and 
|woa — 1^1 1\ 3> Oq^; and (iii) if the detuning of the first 
pulse with the |0) <-> |a) transition and the detuning of 
the second pulse with the |1) <-> |a) transition are the 
same wo a — u>i = u>i a — u>u = A, which creates a two- 
photon Raman resonance between the levels |0) and |1) 
(a standard technique in quantum optics [14,15]) [Fig. 
1(a)], the effective interaction Hamiltonian in the inter- 
action picture with respect to H s = J2i=o 1 a \i) {i\ can 
be written as 

Hi = h{n I Qa e- iAt+i ^ |0) (o| + n[ I a e- tAt+ ^ 1 ' |1) (a\) + h.c. 

(3) 



In above, ui. 
between the levels 



tJ = (Ej — Ei)/h is the transition frequency 
,z)and | j> , and fi* = | ^ <»| * | j> | 
is the Rabi flopping frequencies between the levels |z) 
and \j) (on resonance) generated by the fcth pulse (here, 
i,j = 0,1, a; i^j; 3>* w = J S B* ■ dS; and k = /, II). 

To simplify our presentation, we now replace Qqo and 
Q{ ! a of Eq. (3) by Slj and fljj, respectively. From Eq. 
(3), it is easy to show that when the condition A >> 
fliji is satisfied (i.e., large detuning), the Hamiltonian 
(3) reduces to the following effective Hamiltonian after 
eliminating the intermediate level \a) 



n 2 
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in a rotating frame U — e %ri with rj = — A (coo + en) t, 
where the operators croo = |0) (0| , an = |1) (1| , c~ = 
|0) (1|, and a+ = |1) (0| . In Eq. (4), the first two terms 
are ac-Stark shifts of the levels |0) and |1) , which are 
induced by the two microwave pulses, respectively; while 
the last two terms indicate that the levels |0) and |1) are 
coupled to each other with an effective coupling constant 
g = fli fljj/ A, due to the two-photon Raman transition. 
It is clear that Eq. (4) can be rewritten as 



H, 
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where ujo = (fij — ^j/)/(2A) and a z = an — er o- 

In the following, the two logical states of a SQUID 
qubit are represented by the two lowest levels |0) and |1) 
of the SQUID. The dynamics of the qubit is governed by 
the Hamiltonian (4) or (5). Here, it should be mentioned 



that a term —HA (aoo + fn) m Eq. (4) and another term 

— h n '^" (ooo+o"ii) m Eq. (5) have been safely omitted, 
because these two terms are proportional to the identity 
I = |0) (0| + |1) (1| in the qubit Hilbert space, i.e. they 
bring a common phase factor to the states |0) and |1) 
during the time evolution. 

By means of Eq. (5), one can easily obtain the time 
evolution of the logical states |0) and |1) as follows 



|0) -» (cos ^Jg 2 +oj 2 t + 



too 



: sin + cj 2 t) |0) 
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sin ^Jg 2 + uj 2 t |1) , 



smJg 2 +u 2 t\Q) 



+ (cos Jg 2 + u J 2 t-z ^L_^ sin J g 2 + Z 2 t) |1) . (6) 

v g 2 + 



In particular, for the case of luq = 0, i.e., Qi = flu, or 
ujo <C g, i.e., SQ <C fi/fi/j/fi (here, 5fl = Qj — flu and 
H = Et±0i£) 5 the state rotation is given by 



|0) -» cos gt |0) +ie" 



t>i—<t>ii) 



sin gt |1) , 



|1) 



le 



i(4>l-<pn) 



sin gt |0) + cos gt |1) . 



(7) 



Eq. (6) or (7) shows that the intermediate level \a) is 
not populated during the time evolution. This is due to 
the fact that the level \a) is not involved in the above 
Hamiltonian (4) or (5). 

Based on Eq. (7), one can see that (i) A flipping be- 
tween |0) and |1) (i.e., 7r-rotation) can be implemented 
by setting t — n/(2g); (ii) The Hadamard transforma- 
tion |0) -» (|1) + |0)) /V2 and |1) -» (|1) - |0» /V2 can 
be achieved with t = (tt/2) /(2g) and 4>i — 4>u = tt/2; 
(iii) More importantly, the single-qubit phase shift gate 
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-i(4>i — <fiu) 



|0) and |1) -> e 



i(4>i — <Pn) 



1) can be real- 



ized by a two-step operation: setting the operation time 
t = n/(2g) with an arbitrary (j>j — <f>u for the first step 
and t = 7r/ (2g) and <j)j — <pu = n for the second step, re- 
spectively. Thus, a complete set of single-qubit arbitrary 
operations can be achieved by combining the qubit-state 
rotation (7) with a single-qubit phase shift [13]. 

It is interesting to note that the gate operations de- 
scribed above can be much faster than the conventional 
gate operations via resonant Rabi oscillations between 
the qubit levels |0) and |1) [Fig. 1(a)]. For the sake of 
concreteness, let's choose a SQUID with the following 
parameters: (3l — 1-20, Z — 50fl,ujLc = 5 x 10 11 rad/s 
(i.e., C = 40 £F, L = 100 pH, I c = 3.95 //A), and 
$ x = -0.501*o, where L = 2ttLI c /<P , Z = y/L/C, 
and ujlc = 1/VLC are, respectively, the SQUID's po- 
tential shape parameter, characteristic impedance, and 
characteristic frequency. For simplicity, denote <pij = 
(i\ $ \ j) /$o- A simple numerical calculation shows </>o a = 
8.4xl0" 5 , (f>ia = 5.4xl0~ 5 , and ^oi = 7.9xl0" 7 . Based 
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on the values of these coupling matrix elements, assum- 
ing the amplitudes of the microwave pluses for the above 
three-level gate and the conventional two-level gate on 
the same order, we have 0/ ~ n H ~ 10 2 fi i. Thus, if 
we set the detuning A = 10 max (fij, fi/j) , the effective 
Rabi-flopping frequency g in equation (7) and the Rabi- 
flopping frequency Qqi m the conventional gate operation 
have the relationship g ~ lOfioi, i- e -7 the gate speed in 
the present scheme can be as about ten times as that 
in the conventional scheme. Here, we should point out 
that setting A = 10f2/ or IOQjj is reasonable, because 
if A = 10S1/, the maximum population of the level \a) 
would be about Qj /(Qj + lOOQj) ~ 0.01 or 1% during 
the gate operations. 

Some points may need to be addressed here. Firstly, 
since both the Qi and Qn couplings are far-off resonance, 
no transition between any two levels is induced by a sin- 
gle pulse. Secondly, in order to avoid two-photon Raman 
transitions between the levels |0) and \a) or the levels |1) 
and | a), the conditions u>i ± u>u ^ u>o a , u>i a must be sat- 
isfied. Fortunately, the condition cjj — uu ^ uoa, oj\ a is 
ensured because of uiQ a ,uJi a 3> Woi = wj — u>u. On the 
other hand, the condition uji+ujh ^ ujQ a , ^>ia can be eas- 
ily satisfied by adjusting the SQUID's level structure [see 
Fig. 1(a)]. In fact, all of the conditions required can be 
readily achieved experimentally by adjusting microwave 
frequencies, or level spacings of the SQUID qubit. 

Recently, geometric quantum computing has been paid 
much attention because apart from its fundamental in- 
terest, geometric phase is intrinsically fault-tolerant to 
certain types of computational errors due to its geomet- 
ric property [3,16-20]. On one hand, several basic ideas 
of geometric computing based on Berry phase have been 
proposed by using NMR [16], superconducting electron 
box [17], or trapped ions [18]; and an experimental real- 
ization of the conditional adiabatic phase shift has been 
reported with NMR technique [16]. On the other hand, 
geometric computing based on Aharonov and Anandan 
(A- A) phase [21] is currently getting considerable atten- 
tion and several proposals have been also presented re- 
cently [3,19,20]. In A-A phase geometric manipulation, if 
the external field is perpendicular to the evolution path, 
there is no dynamic phase accumulated in the whole pro- 
cess; thus no extra operation is required to eliminate 
the dynamic phase. In contrast, geometric manipulation 
based on Berry phase always involves the dynamic phase 
and an extra operation is needed to remove its adverse 
effect [3]. In the following, we show explicitly how to per- 
form an arbitrary rotation and a phase shift for a SQUID 
qubit via the A-A phase geometric manipulation [22] . 

We start by rewriting the Hamiltonian (5) as 

H e ff = fkJ Q a z - Kg cos (4>i - fin) a x - fig sin {fii - fin) a y 
= B-a, (8) 

where the fictitious field B= 



{—fig cos (4>i — <fru) , — fig sin (fii — fin) , hu>o} and a are 
the Pauli operators. To see how the logical states |0) and 
|1) are rotated via geometric means, let us first show how 
the eigenstates |±) of a y , defined by a y |±) = ± |±) , will 
undergo a cyclic evolution and obtain a A-A phase under 
the following operations: (i) turn on the two microwave 
pulses / and II, which have fii — fin = (2n + l)n (n is 
an integer). After a time t = n/ {2y/iv 2 + g 2 ), the state 

|+) rotates around the fictitious field B= {fig, 0, Uloq} , 
from |+) in the e y direction to |— ) in the -e y direction 
along curve ABC on the Bloch sphere [see Fig. 1(b)]. (ii) 
Change the two-pulse phase difference to (pi — fin = 2nn. 
After another time t — n/ (2-^/Wq + g 2 ), the state | — ) ro- 
tates back to |+) around B= {— fig, 0, Hloq} along curve 
CDA on the Bloch sphere. In the above operations, the 
trace of the state vector encloses an area on the Bloch 
sphere. After the above cyclic evolution, the state |+) be- 
comes e~ l2t> |+) with a A-A geometric phase —28, where 
29 = 2 arctan(w /5) is the solid angle subtended by the 
area of ABCDA. It should be mentioned that because 
the state vector is always perpendicular to the effec- 
tive magnetic field during the above operations, no dy- 
namical phase is accumulated during the evolution. In 
a similar way, the state |— ) will become e l26 \— ) with 
an accumulated A-A phase of 29 after the above oper- 
ations. Finally, based on |0) = — i (|+) — |— )) j\f2 and 
|1) = (|+) + |— )) /V2, one can see that at the end of the 
above operation, the logical states |0) and |1) are rotated 
as 

|0) — > cosa |0) — sina |1) , 

|1) -» sina|0) +cosa|l) , (9) 

where a = 29 is within the range of [0, 2ir] , i.e., (9) ac- 
complishes a rotation with angle a. 

For the case of f2/ = flu, i.e., uj = 0, the 
Hamiltonian (8) reduces to H e ff =B • & with 

B= { — fig cos ((f>i — fin) ,— fig sin (fii — fin) ,0} . To im- 
plement a single-qubit phase shift gate, one can first 
apply two pulses with a certain phase difference 5fio. 
After a time t — An/ (2^ljiljj) , the logical states 
|0) and |1) rotate around the effective magnetic field 
{— HgcosSfio, — figs'mSfio, 0} to |1) and |0) , respectively. 
In the second step of the operation, one changes the 
phase difference to —Sfio- After another time t = 
An/ (2f2/0//), the logical states rotate around the effec- 
tive magnetic field {—fig cos Sfio , fig sin Sfio , 0} and return 
to the original states |0) and |1) . After this cyclic evolu- 
tion, the logical states acquire geometric A-A phases as 
|0) -» e l<3 |0) and |1) -> |1) , where j3 = n - 28fi a 
is the solid angle subtended by the area of the loop. 
According to [13], combining this phase shift gate with 
the qubit state rotation (9) constitutes a complete set of 
single-qubit arbitrary operations. 

One significant point needs to be made here. For the 
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rotation gate, it is not necessary to change suddenly the 
two-pulse phase difference after the first-step operation. 
Imagine that after the first step, the two pulses are turned 
off for a time interval St in order to adjust the phase 
difference to —54>o- In this case, by a simple calculation, 
we find that after the whole operation, the two logical 
states evolve into 

|0) -> cos a |0) - e-^-^Wsinall) , 

|1) -» sina|0) + e -HEi-Eo)St/h cosa |!) , ( 10 ) 

where the relative phase e - % ( E i- E o)St/h between |0) and 
|1) is induced due to the free evolution of the logical 
states during the period of adjusting the phase difference. 
However, from Eq. (10), one can see that if the time 
interval satisfies St = 2mr/ ( El ^ Ea ) {n is an integer), the 
qubit-state rotation described by Eq. (10) is the same as 
that given by Eq. (9) exactly. 

In summary, we have explicitly shown how to perform 
an arbitrary-operation gate with a SQUID qubit. The 
scheme has three distinct advantages: (i) The intermedi- 
ate level is unpopulated during the gate operation, thus 
decoherence due to energy relaxation from the level \a) 
is minimized; (ii) No tunneling between the qubit lev- 
els |0) and |1) is required, therefore the decay from the 
level |1) can be made negligibly small via increasing the 
potential barrier between the qubit levels; (iii) For the 
3-level gates described here, stronger microwave field can 
be used to achieve faster operation than the conventional 
2-level gates. More interestingly, as shown above, an ar- 
bitrary quantum logic via global A-A phase geometric 
means can be achieved in a SQUID qubit, by extending 
the present method. In addition, the method can also 
be applied to any other type of qubits with a A-type 
level configuration. The present proposal provides new 
approaches to demonstrate a single-qubit arbitrary gate 
with superconducting devices, and we hope that the pro- 
posed approaches will stimulate further theoretical and 
experimental activities in this area. 

Before we conclude, it should be noticed that recently, 
M. H. S. Amin et al. [23] have proposed a similar scheme 
for obtaining a complete set of one-qubit gate with a 
SQUID qubit. In their proposal, an arbitrary gate oper- 
ation is achieved based on a kind of rotation that is per- 
formed over a half period of Rabi oscillation by applying 
microwave to induce transition to the intermediate level 
(a higher state). As addressed in Ref. [23], the proba- 
bility of finding the system in the intermediate level is 
zero again after a half of Rabi period. However, under 
their assumption of small detuning, population in the in- 
termediate level is finite during the evolution, resulting 
in higher probability of spontaneous decay. Finally, we 
should point out that qubit operations via adiabatic pas- 
sage [24] in general are rather slow and require precise 
control over the amplitude of external fields. 
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Figure Captions 

FIG. 1. (a) Simplified level diagram of an rf SQUID 
with A-type three levels |0),|1) and |o). (b) The path 
evolution of the SQUID qubit initially in the state |+) . 
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